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Abstract 

Logical hidden Markov models (LOHMMs) upgrade traditional hidden Markov models 
to deal with sequences of structured symbols in the form of logical atoms, rather than flat 
characters. 

This note formally introduces LOHMMs and presents solutions to the three central in- 
ference problems for LOHMMs: evaluation, most likely hidden state sequence and param- 
eter estimation. The resulting representation and algorithms are experimentally evaluated 
on problems from the domain of bioinformatics. 

1. Introduction 

Hidden Markov models (HMMs) (Rabiner & Juang, 1986) are extremely popular for an- 
alyzing sequential data. Application areas include computational biology, user modelling, 
speech recognition, empirical natural language processing, and robotics. Despite their suc- 
cesses, HMMs have a major weakness: they handle only sequences of flat, i.e., unstruc- 
tured symbols. Yet, in many applications the symbols occurring in sequences are struc- 
tured. Consider, e.g., sequences of UNIX commands, which may have parameters such 
as emacs lohmms.tex, Is, latex lohnuns.tex, . . .Thus, commands are essentially structured. 
Tasks that have been considered for UNIX command sequences include the prediction of 
the next command in the sequence (Davison & Hirsh, 1998), the classification of a command 
sequence in a user category (Korvemaker & Greiner, 2000; Jacobs & Blocked, 2001), and 
anomaly detection (Lane, 1999). Traditional HMMs cannot easily deal with this type of 
structured sequences. Indeed, applying HMMs requires either 1) ignoring the structure of 
the commands (i.e., the parameters), or 2) taking all possible parameters explicitly into 
account. The former approach results in a serious information loss; the latter leads to a 
combinatorial explosion in the number of symbols and parameters of the HMM and as a 
consequence inhibits generalization. 

The above sketched problem with HMMs is akin to the problem of dealing with struc- 
tured examples in traditional machine learning algorithms as studied in the fields of in- 
ductive logic programming (Muggleton & De Raedt, 1994) and multi-relational learn- 
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ing (Dzeroski &: Lavrac, 2001). In this paper, we propose an (inductive) logic programming 
framework, Logical HMMs (LOHMMs), that upgrades HMMs to deal with structure. The 
key idea underlying LOHMMs is to employ logical atoms as structured (output and state) 
symbols. Using logical atoms, the above UNIX command sequence can be represented 
as emacs(lohiiuiis.tex). Is, latex(lohmnis.tex), . . . There are two important motivations for 
using logical atoms at the symbol level. First, variables in the atoms allow one to make 
abstraction of specific symbols. E.g., the logical atom eiiiacs(X, tex) represents all files X 
that a LTgX user tex could edit using emacs. Second, unification allows one to share in- 
formation among states. E.g., the sequence emacs(X, tex), latex(X, tex) denotes that the 
same file is used as an argument for both Emacs and LTf|]X. 

The paper is organized as follows. After reviewing the logical preliminaries, we introduce 
LOHMMs and define their semantics in Section 3; in Section 4, we upgrade the basic 
HMM inference algorithms for use in LOHMMs; we investigate the benefits of LOHMMs in 
Section 5: we show that LOHMMs are strictly more expressive than HMMs, that they can 
be — by design — an order of magnitude smaller than their corresponding propositional 
instantiations, and that unification can yield models, which better fit the data. In Section 6, 
we empirically investigate the benefits of LOHMMs on real world data. Before concluding, 
we discuss related work in Section 7. Proofs of all theorems can be found in the Appendix. 

2. Logical Preliminaries 

A first-order alphabet S is a set of relation symbols r with arity m > 0, written r/m, and a 
set of functor symbols f with arity n > 0, written f/n. If n = then f is called a constant, 
if m = then p is called a propositional variable. (We assume that at least one constant 
is given.) An atom r(t]^, . . . ,t„) is a relation symbol r followed by a bracketed n-tuple of 
terms tj. A term t is a variable V or a functor symbol f (t^, . . . , t;.) immediately followed by 
a bracketed fc-tuple of terms tj. Variables will be written in upper-case, and constant, func- 
tor and predicate symbols lower-case. The symbol _ will denote anonymous variables which 
are read and treated as distinct, new variables each time they are encountered. An iterative 
clause is a formula of the form H B where H (called head) and B (called body) are logical 
atoms. A substitution 9 = {Vi/ti, . . . , V„/t„}, e.g. {X/tex}, is an assignment of terms tj 
to variables Vj. Applying a substitution cr to a term, atom or clause e yields the instanti- 
ated term, atom, or clause ea where all occurrences of the variables Vj are simultaneously 
replaced by the term tj, e.g. ls(X) emacs(F, X){X/tex} yields Is(tex) ^ emacs(F, tex). 
A substitution a is called a unifier for a finite set S of atoms if Sa is singleton. A unifier 9 
for S is called a most general unifier (MGU) for S if, for each unifier a of S, there exists a 
substitution 7 such that a = 9'y. A term, atom or clause E is called ground when it contains 
no variables, i.e., vars{E) = 0. The Herbrand base of S, denoted as hbs, is the set of all 
ground atoms constructed with the predicate and functor symbols in S. The set Ge(A) of 
an atom A consists of all ground atoms A9 that belong to hbs. 

3. Logical Hidden Markov Models 

The logical component of a traditional HMM corresponds to a Mealy machine (Hopcroft 
& Ullman, 1979), i.e., a finite state machine where the output symbols are associated with 
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transitions. This is essentially a propositional representation because the symbols used to 
represent states and output symbols are flat, i.e. not structured. The key idea underlying 
LOHMMs is to replace these flat symbols by abstract symbols. An abstract symbol A is — 
by definition — a logical atom. It is abstract in that it represents the set of all ground, i.e., 
variable- free atoms of A over the alphabet S, denoted by Gx;(A). Ground atoms then play 
the role of the traditional symbols used in a HMMs. 

Example 1 Consider the alphabet Si which has as constant symbols tex, dvi, lunml, 
and lohininl, and as relation symbols emacs/2, ls/1, xdvi/1, latex/2. Then the atom 
emacs(File, tex) represents the set {emacs(hininl, tex), emacs(lohininl, tex)}. We assume 
that the alphabet is typed to avoid useless instantiations such as eniacs(tex, tex)j. 

The use of atoms instead of flat symbols allows us to analyze logical and structured sequences 
such as eniacs(hmml, tex), latex(hniiiil, tex), xdvi(hmml, dvi). 

Definition 1 Abstract transition are expressions of the form p : K B where p € [0, 1], 
and E, B and are atoms. All variables are implicitly assumed to be universally quantified, 
i.e., the scope of variables is a single abstract transition. 

The atoms H and B represent abstract states and represents an abstract output symbol. 
The semantics of an abstract transition p : H <^ B is that if one is in one of the states in 
Gi;(B), say B^B) one will go with probability p to one of the states in Gs(H0b)) say H^b^H; 
while emitting a symbol in Ge(D^b^h)) say O^b^h^o- 

Example 2 Consider c = 0.8 : xdvi(File, dvi) < ^- ^ latex(File, tex). In general 

H, B and do not have to share the same predicate. This is only due to the na- 
ture of our running example. Assume now that we are in state latex(hmnil, tex), i.e. 
9-B, = {File/hmml}. Then c specifies that there is a probability of 0.8 that the next state 
will be in Gs^ (xdvi(hmnil, dvi)) = {xdvi(hmnil, dvi)} ( i.e., the probability is 0.8 that the 
next state will be xdvi(hmnil, dvi)/, and that one of the symbols in Gs^ (latex(hmnil)) = 
{latex(hmml)} ( i.e., latex(hmml)j will be emitted. Abstract states might also be more 
complex such as latex(f ile(FileSteiii, FileExtension), User) 

The above example was simple because and were both empty. The situation be- 
comes more complicated when these substitutions are not empty. Then, the resulting 
state and output symbol sets are not necessarily singletons. Indeed, for the transi- 

, , , latex(File) , , 

tion 0.8 : emacs(File , dvi) < latex(File, tex) the resultmg state set would be 

Gs^ (emacs(File', dvi)) = {einacs(hmnil, tex), emacs(lohmnil, tex)}. Thus the transition 
is non-deterministic because there are two possible resulting states. We therefore need a 
mechanism to assign probabilities to these possible alternatives. 

Definition 2 The selection distribution ^ specifies for each abstract state and observation 
symbol A over the alphabet S a distribution | A) over Gj:,{A). 

To continue our example, let /x(eiiiacs(hmnil, tex) | eniacs(File', tex)) = 0.4 and 
//(emacs(lohinml, tex) | eiiiacs(File', tex)) = 0.6. Then there would be a probabil- 
ity of 0.4 X 0.8 = 0.32 that the next state is emacs(hiiuiil, tex) and of 0.48 that it is 
emacs(lohmnil, tex). 
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Taking /i into account, the meaning of an abstract transition p : E i — B can be sum- 
marized as fohows. Let BOb € Gs(B), EOBOn € GsiEOB) and 09^6^00 G Gj:{O0b0e)- Then the 
model makes a transition from state B^b to H^b^h and emits symbol O^b^h^o with probability 

p ■ I H^b) • ^(O^B^H^O I O^b^h). (1) 

To represent /i, any probabilistic representation can - in principle - be used, e.g. a Bayesian 
network or a Markov chain. Throughout the remainder of the present paper, however, 
we will use a naive Bayes approach. More precisely, we associate to each argument of a 
relation r/m a finite domain D^''™ of constants and a probability distribution p^^^ over 

D^^"^. Let vars(A) = {Vi, . . . , V/} be the variables occurring in an atom A over r/m, and 
let a = {Vi/si, . . . V;/s/} be a substitution grounding A. Each Vj is then considered a 

random variable over the domain -D^^^y ) of the argument arg(Vj) it appears first in. Then, 
fi{Aa I A) = nj=i -^arg(v -^'S- Ai(eniacs(hininl, tex) | emacs(F,E)), is computed as the 

product of P™^'^^^^(hmml) and P2'^^'^^^'^ (tex). 

Thus far the semantics of a single abstract transition has been defined. A LOHMM 
usually consists of multiple abstract transitions and this creates a further complication. 

, , emacs(Flle) , , 

Example 3 Consider 0.8 : latex(File, tex) i emacs(File, tex) ana 

. emacs(Flle) , . 

0.4 : dvi(File) i emacs(File,User). These two abstract transitions make 

conflicting statements about the state resulting from emacs(hmml, tex). Indeed, according 
to the first transition, the probability is 0.8 that the resulting state is latex(hmml, tex) and 
according to the second one it assigns 0.4 to xdvi(hmml). 

There are essentially two ways to deal with this situation. On the one hand, one might want 
to combine and normalize the two transitions and assign a probability of | respectively | . 
On the other hand, one might want to have only one rule firing. In this paper, we chose the 
latter option because it allows us to consider transitions more independently, it simplifies 
learning, and it yields locally interpretable models. We employ the subsumption (or gen- 
erality) relation among the B-parts of the two abstract transitions. Indeed, the B-part of 
the first transition Bi = emacs(File, tex) is more specific than that of the second transi- 
tion B2 = emacs(File, User) because there exists a substitution 6 = {User/tex} such that 
B2O = Bi, i.e., B2 subsumes Bi. Therefore Gs^(Bi) C Ge^(B2) and the first transition can 
be regarded as more informative than the second one. It should therefore be preferred over 
the second one when starting from emacs(hmnil, tex). We will also say that the first tran- 
sition is more specific than the second one. Remark that this generality relation imposes a 
partial order on the set of all transitions. These considerations lead to the strategy of only 
considering the maximally specific transitions that apply to a state in order to determine 
the successor states. This implements a kind of exception handling or default reasoning 
and is akin to Katz's (1987) back-off n-gram models. In back-off n-gram models, the most 
detailed model that is deemed to provide sufficiently reliable information about the current 
context is used. That is, if one encounters an n-gram that is not sufficiently reliable, then 
back-off to use an (n — l)-gram; if that is not reliable either then back-off to level n — 2, etc. 

The conflict resolution strategy will work properly provided that the bodies of all max- 
imally specific transitions (matching a given state) represent the same abstract state. This 
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latex(F) : 0.6 



Figure 1: A logical hidden Markov model. 



can be enforced by requiring the generality relation over the B-parts to be closed under the 
greatest lower bound (gib) for each predicate, i.e., for each pair Bi,B2 of bodies, such that 
6 = mgu(Bi,B2) exists, there is another body B (called lower bound) which subsumes Bi^ 
(therefore also B2O) and is subsumed by Bi,B2, and if there is any other lower bound then 
it is subsumed by B. E.g., if the body of the second abstract transition in our example is 
emacs(hiiiml, User) then the set of abstract transitions would not be closed under gib. 

Finally, in order to specify a prior distribution over states, we assume a finite set T of 
clauses of the form p : E start using a distinguished start symbol such that p is the 
probability of the LOHMM to start in a state of Gs(H). 

By now we are able to formally define logical hidden Markov models. 

Definition 3 A logical hidden Markov model (LOHMM) is a tuple (S,^, A,T) where E is 
a logical alphabet, fj. a selection probability over S, A is a set of abstract transitions, and T 
is a set of abstract transitions encoding a prior distribution. Let B be the set of all atoms 
that occur as body parts of transitions in A. We assume B to be closed under gib and require 

VB G B : y „ p = 1.0 (2) 

and that the probabilities p of clauses in T sum up to 1.0 . 

HMMs are a special cases of LOHMMs in which T, contains only relation symbols of arity 
zero and the selection probability is irrelevant. Thus, LOHMMs directly generalize HMMs. 

LOHMMs can also be represented graphically. Figure 1 contains an example. The under- 
lying language S2 consists of Si together with the constant symbol other which denotes a 
user that does not employ L^T[^. In this graphical notation, nodes represent abstract states 
and black tipped arrows denote abstract transitions. White tipped arrows are used to repre- 
sent meta knowledge. More precisely, white tipped, dashed arrows represent the generality or 
subsumption ordering between abstract states. If we follow a transition to an abstract state 
with an outgoing white tipped, dotted arrow then this dotted arrow will always be followed. 
Dotted arrows are needed because the same abstract state can occur under different cir- 

„ . , , . . / , ,N latex(File) , , 

cumstances. Consider the transition p : latex(File , User J i latex(File, User). 
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start^ ein(F,U) -^em(f i , t)-- -em(F, t) la(F, t) ^^»la(fi,t)^ 

abstract state state abstract state abstract state state 0.6 

0.4 0.7 la(fi) 

Is , , M , em(f2) u. 

^ls{t>^ ls{U') ^ em(f2,o)^ em(F,U) - em(F',U) 

state abstract state state abstract state abstract state 

Figure 2: Generating the observation sequence emacs(hinml), latex(hminl), 
emacs(lohininl), Is by the LOHMM in Figure 1. The command emacs is 
abbreviated by em, f i denotes the filename hmml, f 2 represents lohmml, t denotes 
a tex user, and o some other user. White tipped solid arrows indicate selections. 



em(f 1) 



Even though the atoms in the head and body of the transition are syntactically different they 
represent the same abstract state. To accurately represent the meaning of this transition we 
cannot use a black tipped arrow from latex(File, User) to itself, because this would actu- 

/ N latex(File) , 

ally represent the abstract transition p : latex(File, User) < latex(File, User). 

Furthermore, the graphical representation clarifies that LOHMMs are generative mod- 
els. Let us explain how the model in Figure 1 would generate the observation sequence 
emacs(hniiiil), latex(hmml), emacs(lohmml). Is (cf. Figure 2). It chooses an initial ab- 
stract state, say eniacs(F,U). Since both variables F and U are uninstantiated, the model 
samples the state emacs(hmml, tex) from using ^. As indicated by the dashed ar- 
row, emacs(F, tex) is more specific than emacs(F,U). Moreover, emacs(hinml, tex) matches 
emacs(F, tex). Thus, the model enters emacs (F, tex). Since the value of F was already 
instantiated in the previous abstract state, emacs(hmml, tex) is sampled with probability 
1.0. Now, the model goes over to latex(F, tex), emitting emacs(hmml) because the abstract 
observation emacs (F) is already fully instantiated. Again, since F was already instantiated, 
latex(hmml, tex) is sampled with probability 1.0. Next, we move on to emacs(F',U), emit- 
ting latex(hmml). Variables F' and U in emacs(F',U) were not yet bound; so, values, say 
lohmml and others, are sampled from /x. The dotted arrow brings us back to emacs(F,U). 
Because variables are implicitly universally quantified in abstract transitions, the scope of 
variables is restricted to single abstract transitions. In turn, F is treated as a distinct, 
new variable, and is automatically unified with F', which is bound to lohmml. In contrast, 
variable U is already instantiated. Emitting emacs(lohmml), the model makes a transition 
to Is(U'). Assume that it samples tex for U'. Then, it remains in Is(U') with probability 
0.4 . Considering all possible samples, allows one to prove the following theorem. 

Theorem 1 (Semantics) A logical hidden Markov model over a language T, defines a 
discrete time stochastic process, i.e., a sequence of random variables {Xt)t=i,2,..., where the 
domain of Xt is hb(S) xhb(S). The induced probability measure over the Cartesian product 
(^jhb(S) X hb(S) exists and is unique for each t > and in the limit t — )■ 00. 

Before concluding this section, let us address some design choices underlying LOHMMs. 

First, LOHMMs have been introduced as Mealy machines, i.e., output symbols are 
associated with transitions. Mealy machines fit our logical setting quite intuitively as they 
directly encode the conditional probability -P(0,S'|S) of making a transition from S to S' 
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emitting an observation 0. Logical hidden Markov models define this distribution as 
P(0,S'|S) = V p ■ /i(S' I Hctb) ■ //(O I OVeaH) 

^-^p:^i B 

0' 

where the sum runs over all abstract transitions H i — B such that B is most specific for S. 
Observations correspond to (partially) observed proof steps and, hence, provide information 
shared among heads and bodies of abstract transitions. In contrast, HMMs are usually 
introduced as Moore machines. Here, output symbols are associated with states implicitly 
assuming D and S' to be independent. Thus, P(0,S' | S) factorizes into P(0 | S) • P(S' | S). 
This makes it more difficult to observe information shared among heads and bodies. In 
turn, Moore-LOHMMs are less intuitive and harder to understand. For a more detailed 
discussion of the issue, we refer to Appendix B where we essentially show that - as in the 
propositional case - Mealy- and Moore-LOHMMs are equivalent. 

Second, the naive Bayes approach for the selection distribution reduces the model com- 
plexity at the expense of a lower expressivity: functors are neglected and variables are 
treated independently. Adapting more expressive approaches is an interesting future line of 
research. For instance, Bayesian networks allow one to represent factorial HMMs (Ghahra- 
mani & Jordan, 1997). Factorial HMMs can be viewed as LOHMMs, where the hidden 
states are summarized by a 2 • /c-ary abstract state. The first k arguments encode the k 
state variables, and the last k arguments serve as a memory of the previous joint state, 
of the i-th argument is conditioned on the i + k-th. argument. Markov chains allow one to 
sample compound terms of finite depth such as s(s(s(0))) and to model e.g. misspelled 
filenames. This is akin to generalized HMMs (Kulp, Haussler, Reese, & Eeckman, 1996), in 
which each node may output a finite sequence of symbols rather than a single symbol. 

Finally, LOHMMs - as introduced in the present paper - specify a probability distri- 
bution over all sequences of a given length. Reconsider the LOHMM in Figure 1. Al- 
ready the probabilities of all observation sequences of length 1, i.e.. Is, emacs(hiiuiil), and 
emacs(lohmml)) sum up to 1. More precisely, for each t > it holds that ^ P{Xi = 

xi, . . . , Xt = xt) = l-O ■ In order to model a distribution over sequences of variable length, 
i.e., J2t>o xt ^i-^i ~ ^1' ■ ■ ■ = xt) = 1.0 we may add a distinguished end state. 

The end state is absorbing in that whenever the model makes a transition into this state, 
it terminates the observation sequence generated. 

4. Three Inference Problems for LOHMMs 

As for HMMs, three inference problems are of interest. Let M be a LOHMM and let 
O = Oi, O2, . . . , Or, T > 0, be a finite sequence of ground observations: 

(1) Evaluation: Determine the probability P(0 | M) that sequence was generated by 

the model M. 

(2) Most likely state sequence: Determine the hidden state sequence S* that has most 

likely produced the observation sequence 0, i.e. S* = argmaxs-P(S | 0,M) . 

(3) Parameter estimation: Given a set = {Oi, . . . , Ofc} of observation sequences, de- 

termine the most likely parameters A* for the abstract transitions and the selection 
distribution of M, i.e. A* = argmax;vP(0 | A) . 
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abstract selection abstract selection abstract selection 
transition transition transition 



ls(U': 



em(F,U 




^ ,c\'{\ em(F',o1 ^ ;< 
em(fl,t) ' em(f2,t}^ 

latex(fl,t)-c» latex(fl,t) ' latex(f2,t) 



5o 



5i 



S2 



Figure 3: Trellis induced by the LOHMM in Figure 1. The sets of reachable states at time 
0, 1, . . . are denoted by 5o, 5i, . . . In contrast with HMMs, there is an additional 
layer where the states are sampled from abstract states. 



We will now address each of these problems in turn by upgrading the existing solutions for 
HMMs. This will be realized by computing a grounded trellis as in Figure 3. The possible 
ground successor states of any given state are computed by first selecting the applicable 
abstract transitions and then applying the selection probabilities (while taking into account 
the substitutions) to ground the resulting states. This two-step factorization is coalesced 
into one step for HMMs. 

To evaluate 0, consider the probability of the partial observation sequence Oi, O2, . . . , 0^ 
and (ground) state S at time t, < t < T, given the model M = (S, /i, A, T) 

at{S) := P{0i,02,...,0t,qt = S\ M) 

where qt = S denotes that the system is in state S at time t. As for HMMs, at{S) can be com- 
puted using a dynamic programming approach. For t = 0, we set ao{S) = P{qo = S | M) , 
i.e., ao(S) is the probability of starting in state S and, for t > 0, we compute at{S) based 
on at-i{S'): 



9: 

10: 

11: 



5*0 := {start} 
for t = 1,2, . . . ,r do 
St = 9 

for each S G St^i do 



^* initialize the set of reachable states*/ 

^* initialize the set of reachable states at clock t */ 



return 



foreach maximally specific p:H^B€AuT s.t. (Tb = mgu(S, B) exists do 




foreach S' = HctbCTh S Ge(H(Tb) s.t. 


Ot_i unifies with OctbCTh do 








if S' St then 
St := St U {S'} 
at(S') := 0.0 










at(S') :=at(S') + at_i(S) -p ■ 


;[i(S' 1 HcJb) • /i(Ot-l 1 OCTBfJH) 
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where we assume for the sake of simphcity = start for each abstract transition p : H ^ 
start E T. Furthermore, the boxed parts specify all the differences to the HMM formula: 
unification and fi are taken into account. 

Clearly, as for HMMs P(0 | M) = YlseSr '^T'(S) holds. The computational complexity 
of this forward procedure is 0(T ■ s ■ (|B| + o ■ g)) = 0{T ■ s^) where s = maxt=i^2,...,r I'S'tl , 
o is the maximal number of outgoing abstract transitions with regard to an abstract state, 
and g is the maximal number of ground instances of an abstract state. In a completely 
analogous manner, one can devise a backward procedure to compute 



/3j(S) = P(0t+i,0t_ 



, Ot I = S, M) 



This will be useful for solving Problem (3). 

Having a forward procedure, it is straightforward to adapt the Viterbi algorithm as a 
solution to Problem (2), i.e., for computing the most likely state sequence. Let (5f(S) 
denote the highest probability along a single path at time t which accounts for the first t 
observations and ends in state S, i.e.. 



5t{S)= max P(So,Si, 
So,Si,...,Si_i 



,St_i,St = S,Oi,...,Ot_i|A/) . 



The procedure for finding the most likely state sequence basically follows the forward pro- 
cedure. Instead of summing over all ground transition probabilities in line 10, we maximize 
over them. More precisely, we proceed as follows: 



10 
11 
12 
13 



initialize the set of reachable states*/ 

/* initialize the set of reachable states at clock t */ 



Sq := {start} 

for t = 1,2, ... ,T do 
St = $ 

foreach S G St-i do 

foreach maximally specific p:H-f^B€AuT s.t. ctb = mgu(S, B) exists do 
foreach S' = HctbcTh £ G's(H(Tb) s.t. Ot-i unifies with 0(TbO"h do 
if S' St then 
St := St U {S'} 
5t{S,S') := 0.0 

6t{S,S'):=6t{S,S') + 6t-i{S)-p • //(S' | Hctb) • /i(Ot_i | Octbcth) 
foreach S' G St do 

6t{S') =maxsg5i_i'5t(S,S') 
i)t{S') = argmaxs6St_.i '0t(S,S') 



Here, 5t{S,S') stores the probability of making a transition from S to S' and (with 
■01 (S*) = start for all states S) keeps track of the state maximizing the probability along 
a single path at time t which accounts for the first t observations and ends in state S'. The 
most likely hidden state sequence S* can now be computed as 



and S* 



arg max (5r+i(S) 

SGSt+1 



t+1. 



for t = T,r - 1,... ,1 . 



One can also consider problem (2) on a more abstract level. Instead of considering all 
contributions of different abstract transitions T to a single ground transition from state S 



433 



Kersting, De Raedt, & Raiko 



to state S' in line 10, one might also consider the most likely abstract transition only. This 
is realized by replacing line 10 in the forward procedure with 

at{S') := max(at(S'),at_i(S) - p- /i(S' | Ho-b) ■ n{Qt-i \ OaBan)) • 

This solves the problem of finding the (2') most likely state and abstract transition 
sequence: 

Determine the sequence of states and abstract transitions GT* = 
So, To, 5i, Ti, S2, . . . , St, Tt, St+1 where there exists substitutions 9i with Si+i ^ 
Si = TiOi that has most likely produced the observation sequence 0, i.e. 
GT* = argmaxGT-P(GT | 0,M) . 

Thus, logical hidden Markov models also pose new types of inference problems. 

For parameter estimation, we have to estimate the maximum likelihood transition 
probabilities and selection distributions. To estimate the former, we upgrade the well-known 
Baum- Welch algorithm (Baum, 1972) for estimating the maximum likelihood parameters 
of HMMs and probabilistic context-free grammars. 

For HMMs, the Baum- Welch algorithm computes the improved estimate p of the tran- 
sition probability of some (ground) transition T = p : H 4^ B by taking the ratio 

H'i BGAUT 

between the expected number ^(T) of times of making the transitions T at any time given 
the model M and an observation sequence 0, and the total number of times a transitions 
is made from B at any time given M and 0. 

Basically the same applies when T is an abstract transition. However, we have to be 
a little bit more careful because we have no direct access to ^(T). Let ^t(gcl,T) be the 

GD 

probability of following the abstract transition T via its ground instance gel = p : GH i GB 

at time t, i.e.. 



at{GB) ■ p ■ f3t+i{GE) 



6(gcl, T) = p{Q\ M) ■ /^(™ I "^b) • K^t-i I O^b^h) , (4) 



where (Tb,c7"h are as in the forward procedure (see above) and P{0 \ M) is the probability 
that the model generated the sequence 0. Again, the boxed terms constitute the main 
difference to the corresponding HMM formula. In order to apply Equation (3) to compute 
improved estimates of probabilities associated with abstract transitions, we set 



T T 



e(T)= j;ei(T) = 5^5^6(gcl,T) 



t=\ t=l 



where the inner sum runs over all ground instances of T. 

This leads to the following re-estimation method, where we assume that the sets Si of 
reachable states are reused from the computations of the a- and /3-values: 
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1 
2 
3 
4 
5 
6 



/* initialization of expected counts */ 



foreach T € A U T do 



I ^(T) := m /* or if not using pseudocounts */ 
/* compute expected counts */ 



for t = 0, 1, . . . ,T do 



foreach S G S't do 



7: 



foreach max. specific T = p:H^ — BsAuT s.t. cr^ = mgu(S, B) exists do 



8 



foreach S' = HctbCTh S GE(HcrB) s.t. S' € St+i A mgu(0(, DctbCTh) exists do 



9: 



C(T) := e(T) + at(S) • p • /3t+i(S')/P(0 | Af)- /x(S' | Hctb) • ^(Ot_i | 0(Jb(Th) 



Here, equation (4) can be found in line 9. In line 3, we set pseudocounts as small sample- 
size regularizers. Other methods to avoid a biased underestimate of probabilities and even 
zero probabilities such as m-estimates (see e.g., Mitchell, 1997) can be easily adapted. 

To estimate the selection probabilities, recall that /i follows a naive Bayes scheme. There- 
fore, the estimated probability for a domain element d £ D for some domain D is the ratio 
between the number of times d is selected and the number of times any d' € D is selected. 
The procedure for computing the ^-values can thus be reused. 

Altogether, the Baum- Welch algorithm works as follows: While not converged, (1) es- 
timate the abstract transition probabilities, and (2) the selection probabilities. Since it is 
an instance of the EM algorithm, it increases the likelihood of the data with every update, 
and according to McLachlan and Krishnan (1997), it is guaranteed to reach a stationary 
point. All standard techniques to overcome limitations of EM algorithms are applicable. 
The computational complexity (per iteration) is 0{k ■ (a + d)) = 0{k ■ T ■ s"^ + k ■ d) where 
k is the number of sequences, a is the complexity of computing the a- values (see above), 
and d is the sum over the sizes of domains associated to predicates. Recently, Kersting 
and Raiko (2005) combined the Baum- Welch algorithm with structure search for model 
selection of logical hidden Markov models using inductive logic programming (Muggleton 
& De Raedt, 1994) refinement operators. The refinement operators account for different 
abstraction levels which have to be explored. 

5. Advantages of LOHMMs 

In this section, we will investigate the benefits of LOHMMs: (1) LOHMMs are strictly 
more expressive than HMMs, and (2), using abstraction, logical variables and unification 
can be beneficial. More specifically, with (2), we will show that 

(Bl) LOHMMs can be — by design — smaller than their propositional instantiations, and 
(B2) unification can yield better log-likelihood estimates. 

5.1 On the Expressivity of LOHMMs 

Whereas HMMs specify probability distributions over regular languages, LOHMMs specify 
probability distributions over more expressive languages. 
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Theorem 2 For any (consistent) probabilistic context-free grammar (PCFG) G for some 
language C there exists a LOHMM M s.t. Pg{w) = Pm{w) for all w £ C. 

The proof (see Appendix C) makes use of abstract states of unbounded 'depth'. More 
precisely, functors are used to implement a stack. Without functors, LOHMMs cannot 
encode PCFGs and, because the Herbrand base is finite, it can be proven that there always 
exists an equivalent HMM. 

Furthermore, if functors are allowed, LOHMMs are strictly more expressive than PCFGs. 
They can specify probability distributions over some languages that are context-sensitive: 



1.0 


stack(s(0), s(0)) 




Start 






0.8 


stack(s(X),s(X)) 




stack(X,X) 






0.2 


unstack(s(X), s(X)) 


a 
■<r- 


stack(X,X) 






1.0 


unstack(X, Y) 


b 
■<r- 


unstack(s 


:x) 


,Y) 


1.0 


unstack(s(0), Y) 




unstack(s 


:o) 


,s(Y)) 


1.0 


end 


end 

< 


unstack(s 


:o) 


,s(0)) 



The LOHMM defines a distribution over {a"6"c" | n > 0}. 

Finally, the use of logical variables also enables one to deal with identifiers. Identifiers 
are special types of constants that denote objects. Indeed, recall the UNIX command 
sequence emacs lohmms.tex. Is, latex lohmms.tex, . . . from the introduction. The filename 
lohmms.tex is an identifier. Usually, the specific identifiers do not matter but rather the 
fact that the same object occurs multiple times in the sequence. LOHMMs can easily deal 
with identifiers by setting the selection probability /i to a constant for the arguments in 
which identifiers can occur. Unification then takes care of the necessary variable bindings. 

5.2 Benefits of Abstraction through Variables and Unification 

Reconsider the domain of UNIX command sequences. Unix users oftenly reuse a newly cre- 
ated directory in subsequent commands such as inmkdir(vtlOOx), cd(vtlOOx), Is(vtlOOx) . 
Unification should allow us to elegantly employ this information because it allows us to spec- 
ify that, after observing the created directory, the model makes a transition into a state 
where the newly created directory is used: 

pi : cd(Dir, mkdir) inkdir(Dir, com) and p2 '■ cd(_, mkdir) mkdir(Dir, com) 

If the first transition is followed, the cd command will move to the newly created directory; 
if the second transition is followed, it is not specified which directory cd will move to. Thus, 
the LOHMM captures the reuse of created directories as an argument of future commands. 
Moreover, the LOHMM encodes the simplest possible case to show the benefits of unifica- 
tion. At any time, the observation sequence uniquely determines the state sequence, and 
functors are not used. Therefore, we left out the abstract output symbols associated with 
abstract transitions. In total, the LOHMM U, modelling the reuse of directories, consists 
of 542 parameters only but still covers more than 451000 (ground) states, see Appendix D 
for the complete model. The compression in the number of parameters supports (Bl). 

To empirically investigate the benefits of unification, we compare U with the variant 
of U where no variables are shared, i.e., no unification is used such that for instance the 
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first transition above is not allowed, see Appendix D. has 164 parameters less than U . 
We computed the following zero-one win function 

/(O) = 1^ [\ogPu{0)-\ogPN{0)] > 
I otherwise 

leave-one-out cross- validated on Unix shell logs collected by Greenberg (1988). Overall, 
the data consists of 168 users of four groups: computer scientists, nonprogrammers, novices 
and others. About 300000 commands have been logged with an average of 110 sessions 
per user. We present here results for a subset of the data. We considered all computer 
scientist sessions in which at least a single mkdir command appears. These yield 283 logical 
sequences over in total 3286 ground atoms. The LOO win was 81.63%. Other LOO statistics 
are also in favor of U: 





training 


test 


logP(O) 


Piv(O) 


logP(O) 


lof. ^^^^^ 
Pn(0) 


u 


-11361.0 


1795.3 


-42.8 


7.91 


N 


-13157.0 


-50.7 



Thus, although U has 164 parameters more than A, it shows a better generalization per- 
formance. This result supports (B2). A pattern often found in U was ^ 

0.15 : cd(Dir, mkdir) -(— mkdir(Dir, com) and 0.08 : cd(_, mkdir) mkdir(Dir, com) 

favoring changing to the directory just made. This knowledge cannot be captured in N 

0.25: cd(_, mkdir) ^ mkdir(Dir, com). 

The results clearly show that abstraction through variables and unification can be beneficial 
for some applications, i.e., (Bl) and (B2) hold. 



6. Real World Applications 

Our intentions here are to investigate whether LOHMMs can be applied to real world 
domains. More precisely, we will investigate whether benefits (Bl) and (B2) can also be 
exploited in real world application domains. Additionally, we will investigate whether 

(B3) LOHMMs are competitive with ILP algorithms that can also utilize unification and 
abstraction through variables, and 

(B4) LOHMMs can handle tree-structured data similar to PCFGs. 

To this aim, we conducted experiments on two bioinformatics application domains: protein 
fold recognition (Kersting, Raiko, Kramer, &: De Raedt, 2003) and mRNA signal structure 
detection (Horvath, Wrobel, &: Bohnebeck, 2001). Both application domains are multiclass 
problems with five different classes each. 

1. The sum of probabilities is not the same (0.15 + 0.08 = 0.23 7^ 0.25) because of the use of pseudo counts 
and because of the subhminal non-determinism (w.r.t. abstract states) in U , i.e., in case that the first 
transition fires, the second one also fires. 
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6.1 Methodology 

In order to tackle the multiclass problem with LOHMMs, we followed a plug-in estimate 
approach. Let {ci,C2, . . . ,Cfc} be the set of possible classes. Given a finite set of training 
examples {{xi,yi)}^^i C ^ x {ci, C2, . . . , Cn}, one tries to find f : X ^ {ci, C2, . . . , c^} 

/(x) = arg max P{x | M, A*) • P(c) . (5) 

ce{ci,c2,.--,Cfc} 

with low approximation error on the training data as well as on unseen examples. In 
Equation (5), M denotes the model structure which is the same for all classes, A* denotes 
the maximum likelihood parameters of M for class c estimated on the training examples 
with Ui = c only, and -P(c) is the prior class distribution. 

We implemented the Baum- Welch algorithm (with pseudocounts m, see line 3) for maxi- 
mum likelihood parameter estimation using the Prolog system Yap-4.4.4. In all experiments, 
we set m = 1 and let the Baum- Welch algorithm stop if the change in log-likelihood was 
less than 0.1 from one iteration to the next. The experiments were ran on a Pentium-IV 

3.2 GHz Linux machine. 

6.2 Protein Fold Recognition 

Protein fold recognition is concerned with how proteins fold in nature, i.e., their three- 
dimensional structures. This is an important problem as the biological functions of proteins 
depend on the way they fold. A common approach is to use database searches to find pro- 
teins (of known fold) similar to a newly discovered protein (of unknown fold). To facilitate 
protein fold recognition, several expert-based classification schemes of proteins have been 
developed that group the current set of known protein structures according to the similarity 
of their folds. For instance, the structural classification of proteins (Hubbard, Murzin, Bren- 
ner, & Chotia, 1997) (SCOP) database hierarchically organizes proteins according to their 
structures and evolutionary origin. From a machine learning perspective, SCOP induces a 
classification problem: given a protein of unknown fold, assign it to the best matching group 
of the classification scheme. This protein fold classification problem has been investigated 
by Turcotte, Muggleton, and Sternberg (2001) based on the inductive logic programming 
(ILP) system PROGOL and by Kersting et al. (2003) based on LOHMMs. 

The secondary structure of protein domains^ can elegantly be represented as logical se- 
quences. For example, the secondary structure of the Ribosomal protein L4 is represented as 

st(null, 2), he(right, alpha, 6), st(plus, 2), he(right, alpha, 4), st(plus, 2), 
he(right, alpha, 4), st(plus, 3), he(right, alpha, 4), st(plus, l), he(hright, alpha, 6) 

Helices of a certain type, orientation and length he{HelixType, HelixOrientation, Length), 
and strands of a certain orientation and length st{StrandOrientation, Length) are atoms over 
logical predicates. The application of traditional HMMs to such sequences requires one to 
either ignore the structure of helices and strands, which results in a loss of information, or to 
take all possible combinations (of arguments such as orientation and length) into account, 
which leads to a combinatorial explosion in the number of parameters 

2. A domain can be viewed as a sub-section of a protein which appears in a number of distantly related 
proteins and which can fold independently of the rest of the protein. 
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Figure 4: Scheme of a left-to-right LOHMM block model. 

The results reported by Kersting et al. (2003) indicate that LOHMMs are well-suited 
for protein fold classification: the number of parameters of a LOHMM can by an order of 
magnitude be smaller than the number of a corresponding HMM (120 versus approximately 
62000) and the generalization performance, a 74% accuracy, is comparable to Turcotte 
et al.'s (2001) result based on the ILP system Progol, a 75% accuracy. Kersting et al. 
(2003), however, do not cross- validate their results nor investigate - as it is common in 
bioinformatics - the impact of primary sequence similarity on the classification accuracy. For 
instance, the two most commonly requested ASTRAL subsets are the subset of sequences 
with less than 95% identity to each other (95 cut) and with less than 40% identity to each 
other (40 cut). Motivated by this, we conducted the following new experiments. 

The data consists of logical sequences of the secondary structure of protein domains. As 
in the work of Kersting et al. (2003), the task is to predict one of the five most populated 
SCOP folds of alpha and beta proteins (a/b): TIM beta/alpha-barrel (fold 1), NAD(P)- 
binding Rossmann-fold domains (fold 2), Ribosomal protein L4 (fold 23), Cysteine hydrolase 
(fold 37), and Phosphotyrosine protein phosphatases I-like (fold 55). The class of a/b 
proteins consists of proteins with mainly parallel beta sheets (beta-alpha-beta units). The 
data have been extracted automatically from the ASTRAL dataset version 1.65 (Chandonia, 
Hon, Walker, Lo Conte, P.Koehl, &: Brenner, 2004) for the 95 cut and for the 40 cut. As 
in the work of Kersting et al. (2003), we consider strands and helices only, i.e., coils and 
isolated strands are discarded. For the 95 cut, this yields 816 logical sequences consisting 
of in total 22210 ground atoms. The number of sequences in the classes are listed as 293, 
151, 87, 195, and 90. For the 40 cut, this yields 523 logical sequences consisting of in total 
14986 ground atoms. The number of sequences in the classes are listed as 182, 100, 66, 122, 
and 53. 



LOHMM structure: The used LOHMM structure follows a left-to-right block topology, 
see Figure 4, to model blocks of consecutive helices (resp. strands). Being in a Block of 
some size s, say 3, the model will remain in the same block for s = 3 time steps. A similar 
idea has been used to model haplotypes (Koivisto, Perola, Varilo, Hennah, Ekelund, Lukk, 
Peltonen, Ukkonen, &: Mannila, 2002; Koivisto, Kivioja, Mannila, Rastas, &: Ukkonen, 
2004). In contrast to common HMM block models (Won, Priigel-Bennett, &: Krogh, 2004), 
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the transition parameters are shared within each block and one can ensure that the model 
makes a transition to the next state s{Block) only at the end of a block; in our example 
after exactly 3 intra-block transitions. Furthermore, there are specific abstract transitions 
for all helix types and strand orientations to model the priori distribution, the intra- and 
the inter-block transitions. The number of blocks and their sizes were chosen according 
to the empirical distribution over sequence lengths in the data so that the beginning and 
the ending of protein domains was likely captured in detail. This yield the following block 
structure 

□□ ••• n i i n ••• n i ii i i 

1 2 19 20 27 28 40 41 46 47 61 62 76 77 

where the numbers denote the positions within protein domains. Furthermore, note that 
the last block gathers all remaining transitions. The blocks themselves are modelled using 
hidden abstract states over 

hc{HelixType, HelixOrientation, Length, Block) and sc{StrandOrientation, Length, Block) . 

Here, Length denotes the number of consecutive bases the structure element consists of. 
The length was discretized into 10 bins such that the original lengths were uniformally 
distributed. In total, the LOHMM has 295 parameters. The corresponding HMM without 
parameter sharing has more than 65200 parameters. This clearly confirms (Bl). 

Results: We performed a 10- fold cross-validation. On the 95 cut dataset, the accuracy was 
76% and took approx. 25 minutes per cross-validation iteration; on the 40 cut, the accuracy 
was 73% and took approx. 12 minutes per cross-validation iteration. The results validate 
Kersting et al.'s (2003) results and, in turn, clearly show that (B3) holds. Moreover, the 
novel results on the 40 cut dataset indicate that the similarities detected by the LOHMMs 
between the protein domain structures were not accompanied by high sequence similarity. 

6.3 mRNA Signal Structure Detection 

mRNA sequences consist of bases (guanine, adenine, uracil, cytosine) and fold intramolec- 
ularly to form a number of short base-paired stems (Durbin, Eddy, Krogh, & Mitchison, 
1998). This base-paired structure is called the secondary structure, cf. Figures 5 and 6. The 
secondary structure contains special subsequences called signal structures that are responsi- 
ble for special biological functions, such as RNA-protein interactions and cellular transport. 
The function of each signal structure class is based on the common characteristic binding 
site of all class elements. The elements are not necessarily identical but very similar. They 
can vary in topology (tree structure), in size (number of constituting bases), and in base 
sequence. 

The goal of our experiments was to recognize instances of signal structures classes in 
mRNA molecules. The first application of relational learning to recognize the signal struc- 
ture class of mRNA molecules was described in the works of Bohnebeck, Horvath, and 
Wrobel (1998) and of Horvath et al. (2001), where the relational instance-based learner 
RIBL was applied. The dataset ^ we used was similar to the one described by Horvath 

3. The dataset is not the same as described in the work by Horvath et al. (2001) because we could not obtain 
the original dataset. We will compare to the smaller data set used by Horvath et al., which consisted of 
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et al. (2001). It consisted of 93 mRNA secondary structure sequences. More precisely, it was 
composed of 15 and 5 SECIS (Selenocysteine Insertion Sequence), 27 IRE (Iron Responsive 
Element), 36 TAR (Trans Activating Region) and 10 histone stem loops constituting five 
classes. 

The secondary structure is composed of different building blocks such as stacking region, 
hairpin loops, interior loops etc. In contrast to the secondary structure of proteins that forms 
chains, the secondary structure of mRNA forms a tree. As trees can not easily be handled 
using HMMs, mRNA secondary structure data is more challenging than that of proteins. 
Moreover, Horvath et al. (2001) report that making the tree structure available to RIBL 
as background knowledge had an influence on the classification accuracy. More precisely, 
using a simple chain representation RIBL achieved a 77.2% leave-one-out cross-validation 
(LOO) accuracy whereas using the tree structure as background knowledge RIBL achieved 
a 95.4% LOO accuracy. 

We followed Horvath et al.'s experimental setup, that is, we adapted their data repre- 
sentations to LOHMMs and compared a chain model with a tree model. 



Chain Representation: In the chain representation (see also Figure 5), 
signal structures are described by single{TypeSingle, Position, Acid) or 
helical^TypeHelical, Position, Acid, Acid). Depending on its type, a structure el- 
ement is represented by either single/3 or helical/4. Their first argument 
TypeSingle (resp. TypeHelical) specifies the type of the structure element, i.e., 
single, bulges, bulges, hairpin (resp. stem). The argument Position is the posi- 
tion of the sequence element within the corresponding structure element counted down, 
i.e.^, {n^'^(O), n^^(O), . . . , n^(0)}. The maximal position was set to 13 as this was the 
maximal position observed in the data. The last argument encodes the observed nucleotide 
(pair). 

The used LOHMM structure follows again the left-to-right block structure shown in 
Figure 4. Its underlying idea is to model blocks of consecutive helical structure ele- 
ments. The hidden states are modelled using single(TypeSingle, Position, Acid, Block) 
and helical^TypeHelical, Position, Acid, Acid, Block). Being in a Block of consecutive he- 
lical (resp. single) structure elements, the model will remain in the Block or transition to a 
single element. The transition to a single (resp. helical) element only occurs at Position 
n(0). At all other positions n{Position), there were transitions from helical (resp. single) 
structure elements to helical (resp. single) structure elements at Position capturing the dy- 
namics of the nucleotide pairs (resp. nucleotides) within structure elements. For instance. 



66 signal structures and is very close to our data set. On a larger data set (with 400 structures) Horvath 
et al. report an error rate of 3.8% . 
4. n'"(0) is shorthand for the recursive application of the functor n on m times, i.e., for position m. 
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helical(steni, n(0), c, g 
helical(stem, n(n(0)), c 
helical{stem, n(n(n(0))), c, g 

single{bulge5, n(0), a 
single (bulges, n(n(0)) 
single(bulge5, n(n(n(0}}}, g 

helical(stera, n(0), c, g 
helical{stem, n(n(0)), c, g 

single(bulge5, n(0), a 

helical(stein, n(0) , a, a^ 
helical (stem, n(n(0) ) , u, a 
helical(stem, n(n(n(0))), u, g 
helical(steni, n(n(n(n(0)))) , u, 
helical(steni, n(n(n(n(n(0))))) , c, 
helical(stem, n(n(n(n(n(n(0)))))) , u, 
helical(stem, n(n(n(n(n(n(n(0))))))) , a, 



single(liairpin, n(n(n(0))), a). 
single(hairpin, n(n(0)), u). 
single(hairpin, n(0), u) . 



single (bulges, n(0) , a). 



Figure 5: The chain representation of a SECIS signal structure. The ground atoms are 
ordered clockwise starting with helical(steni, ii(ii(ii(n(n(n(n(0))))))), a, u) at the 
lower left-hand side corner. 



the transitions for block n(0) at posit 

a : he(steni, n(0), X, Y, n(0)) 

b: he(stem,n(0),Y,X,n(0)) 

c: he(stem, n(0), X, _, n(0)) 

d : he(stem, n(0),_, Y, n(0)) 

e: he(stem, n(0), _, _, n(0)) 



on n(n(0)) were 

Pa:he(stein,n(0),X,Y) 



Pi,:he(stem,n(0),X,Y) 



Pc:lie(stem,n(0),X,Y) 



Pd:he(stem,n(0),X,Y) 



Pe:h6(stem,n(0),X,Y) 



he(stem, n(n(0)), X, Y, n(0))) 
he(stem, n(n(0)), X, Y, n(0))) 
he(stem, n(n(0)), X, Y, n(0))) 
he(stem, ii(n(0)), X, Y, ii(0))) 
he(stem, ii(n(0)), X, Y, ii(0))) 

In total, there were 5 possible blocks as this was the maximal number of blocks of consecutive 
helical structure elements observed in the data. Overall, the LOHMM has 702 parameters. 
In contrast, the corresponding HMM has more than 16600 transitions validating (Bl). 

Results: The LOO test log-likelihood was —63.7, and an EM iteration took on average 
26 seconds. 

Without the unification-based transitions b-d, i.e., using only the abstract transitions 

, , , , pa:lie(stem,n(0),X,Y) , / / \n / 

a: he(stem,n(0),X,Y,n(0)) 4- ^ he(stein, ii(n(0)), X, Y, n(0))) 



e : 



he(stem, n(0), _, _, n(0)) 



Pe:lie(stem,n(0),X,Y) 



he(stem, ii(n(0)), X, Y, n(0))), 



the model has 506 parameters. The LOO test log-likelihood was —64.21, and an EM iter- 
ation took on average 20 seconds. The difference in LOO test log-likelihood is statistically 
significant (paired t-test, p = 0.01). 

Omitting even transition a, the LOO test log-likelihood dropped to —66.06, and the 
average time per EM iteration was 18 seconds. The model has 341 parameters. The 
difference in average LOO log-likelihood is statistically significant (paired t-test, p = 0.001). 

The results clearly show that unification can yield better LOO test log-likelihoods, i.e., 
(B2) holds. 
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nucleotide_pair((c, g)) 
nucleotide_pair((c, g)) 
nucleotide_pair((c, g)). 
helical(s(s(s(s(s(0))))), s{s(s{0))), [c] , stem, n{n{n{0)))) 

nucleotide(a) 
nucleotide(a), 
nucleotide(g) . 

single{s{s{s{s{0)))), s{s{s{0))), [],bulge5, n{n{n{0)))) 
nucleotide_pair((c, g)) 
nucleotide_pair((c, g)). 
helical(s(s(s(0))), s{0), [c, c, c] , stem, n(n(0))) 
nucleotide(a) . 
single(s{s(0)), s{0), [] , bulge5, n(0)) 
nucleotide_pair( (a, a)) 
nucleotide_pair((u, a)) 
nucleotide_pair((u, g)) 
nucleotide_pair{{u, a)) 
nucleotide_pair((c, a)) 
nucleotide_pair{{u, a)) 
nucleotide_pair( (a, u)) 
helical{s{0), 0, [c, c] , stem, n(n{n(n{n(n{n(0)))))))) 
root(0, root, [c]) 



single{s{s{s{s{s{s(0)))))), s(s(s(s(s(0))))) , 
[], hairpin, n(n(n(0)))). 

nucleotide{a) . 
nucleotide{u) . 
nucleotide{u) . 



single(s(s(s(s(s(s(s(0))))))), s(s(s(0))), 
[], bulges, n(0)). 

nucleotide(a) . 







s(0) 
/ \ 
s(s(0))s(s(s(0))) 

s(s{s(s{0))))^ I ^s(s(s(s(s(s(s(0))))))) 
s(s(s(s(s(0)))) 



a(s(s(s(s(s(0)))))) 



Figure 6: The tree representation of a SECIS signal structure, (a) The logical sequence, 
i.e., the sequence of ground atoms representing the SECIS signal structure. The 
ground atoms are ordered clockwise starting with root(0,root, [c]) in the lower 
left-hand side corner, (b) The tree formed by the secondary structure elements. 



Tree Representation: In the tree representation (see Figure 6 (a)), the idea is to capture 
the tree structure formed by the secondary structure elements, see Figure 6 (b). Each 
training instance is described as a sequence of ground facts over 

root(0, root, ^Children), 

helical(ID, ParentID, ^Children, Type, Size), 

nucleotide jpa.ir{BasePair) , 

single(ZD, ParentID, ^Children, Type, Size), 

nucleotide(5ase) . 

Here, ID and ParentID are natural numbers 0, s(0), s(s(0)), . . . encoding the child- 
parent relation, ^Children denotes the number^ of children [], [c], [c, c], . . ., Type is the 
type of the structure element such as stem, hairpin, .. ., and Size is a natural number 
0, n(0), n(n(0)), . . . Atoms root{0, root, ^Children) are used to root the topology. The 
maximal ^Children was 9 and the maximal Size was 13 as this was the maximal value 
observed in the data. 

As trees can not easily be handled using HMMs, we used a LOHMM which basically 
encodes a PCFG. Due to Theorem 2, this is possible. The used LOHMM structure can be 
found in Appendix E. It processes the mRNA trees in in-order. Unification is only used for 
parsing the tree. As for the chain representation, we used a Position argument in the hidden 
states to encode the dynamics of nucleotides (nucleotide pairs) within secondary structure 

5. Here, we use the Prolog short hand notation [■] for hsts. A hst either is the constant [] representing the 
empty list, or is a compound term with functor ./2 and two arguments, which are respectively the head 
and tail of the list. Thus [a, b, c] is the compound term .(a, .(b, .(c, []))). 
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elements. The maximal Position was again 13. In contrast to the chain representation, 
nucleotide pairs such as (a, u) are treated as constants. Thus, the argument BasePair 
consists of 16 elements. 

Results: The LOO test log-likelihood was —55.56. Thus, exploiting the tree structure 
yields better probabilistic models. On average, an EM iteration took 14 seconds. Overall, 
the result shows that (B4) holds. 

Although the Baum- Welch algorithm attempts to maximize a different objective func- 
tion, namely the likelihood of the data, it is interesting to compare LOHMMs and RIBL in 
terms of classification accuracy. 

Classification Accuracy: On the chain representation, the LOO accuracies of all 
LOHMMs were 99% (92/93). This is a considerable improvement on RIBL's 77.2% (51/66) 
LOO accuracy for this representation. On the tree representation, the LOHMM also 
achieved a LOO accuracy of 99% (92/93). This is comparable to RIBL's LOO accuracy of 
97% (64/66) on this kind of representation. 

Thus, already the chain LOHMMs show marked increases in LOO accuracy when com- 
pared to RIBL (Horvath et al., 2001). In order to achieve similar LOO accuracies, Horvath 
et al. (2001) had to make the tree structure available to RIBL as background knowledge. 
For LOHMMs, this had a significant influence on the LOO test log- likelihood, but not on 
the LOO accuracies. This clearly supports (B3). Moreover, according to Horvath et al., 
the mRNA application can also be considered a success in terms of the application domain, 
although this was not the primary goal of our experiments. There exist also alternative 
parameter estimation techniques and other models, such as covariance models (Eddy &: 
Durbin, 1994) or pair hidden Markov models (Sakakibara, 2003), that might have been 
used as well as a basis for comparison. However, as LOHMMs employ (inductive) logic pro- 
gramming principles, it is appropriate to compare with other systems within this paradigm 
such as RIBL. 

7. Related Work 

LOHMMs combine two different research directions. On the one hand, they are related to 
several extensions of HMMs and probabilistic grammars. On the other hand, they are also 
related to the recent interest in combining inductive logic programming principles with 
probability theory (De Raedt & Kersting, 2003, 2004). 

In the first type of approaches, the underlying idea is to upgrade HMMs and probabilistic 
grammars to represent more structured state spaces. 

Hierarchical HMMs (Fine, Singer, & Tishby, 1998), factorial HMMs (Ghahramani Sz 
Jordan, 1997), and HMMs based on tree automata (Frasconi, Soda, &: Vullo, 2002) decom- 
pose the state variables into smaller units. In hierarchical HMMs states themselves can be 
HMMs, in factorial HMMs they can be factored into k state variables which depend on one 
another only through the observation, and in tree based HMMs the represented probability 
distributions are defined over tree structures. The key difference with LOHMMs is that 
these approaches do not employ the logical concept of unification. Unification is essential 
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because it allows us to introduce abstract transitions, which do not consist of more detailed 
states. As our experimental evidence shows, sharing information among abstract states by 
means of unification can lead to more accurate model estimation. The same holds for re- 
lational Markov models (RMMs) (Anderson, Domingos, & Weld, 2002) to which LOHMMs 
are most closely related. In RMMs, states can be of different types, with each type described 
by a different set of variables. The domain of each variable can be hierarchically structured. 
The main differences between LOHMMs and RMMs are that RMMs do not either support 
variable binding nor unification nor hidden states. 

The equivalent of HMMs for context-free languages are probabilistic context-free gram- 
mars (PCFGs). Like HMMs, they do not consider sequences of logical atoms and do not 
employ unification. Nevertheless, there is a formal resemblance between the Baum- Welch 
algorithms for LOHMMs and for PCFGs. In case that a LOHMM encodes a PCFG both 
algorithms are identical from a theoretical point of view. They re-estimate the parameters 
as the ratio of the expected number of times a transition (resp. production) is used and the 
expected number of times a transition (resp. production) might have been used. The proof 
of Theorem 2 assumes that the PCFG is given in Greibach normal form^ (GNF) and uses a 
pushdown automaton to parse sentences. For grammars in GNF, pushdown automata are 
common for parsing. In contrast, the actual computations of the Baum- Welch algorithm 
for PCFGs, the so called Inside-Outside algorithm (Baker, 1979; Lari &: Young, 1990), is 
usually formulated for grammars in Chomsky normal formJ . The Inside-Outside algorithm 
can make use of the efficient CYK algorithm (Hopcroft & Ullman, 1979) for parsing strings. 

An alternative to learning PCFGs from strings only is to learn from more structured data 
such as skeletons^ which are derivation trees with the nonterminal nodes removed (Levy & 
Joshi, 1978). Skeletons are exactly the set of trees accepted by skeletal tree automata (STA). 
Informally, an STA, when given a tree as input, processes the tree bottom up, assigning a 
state to each node based on the states of that node's children. The STA accepts a tree iff 
it assigns a final state to the root of the tree. Due to this automata-based characterization 
of the skeletons of derivation trees, the learning problem of (P)CFGs can be reduced to 
the problem of an STA. In particular, STA techniques have been adapted to learning tree 
grammars and (P)CFGs (Sakakibara, 1992; Sakakibara et al., 1994) efficiently. 

PCFGs have been extended in several ways. Most closely related to LOHMMs are 
unification-based grammars which have been extensively studied in computational linguis- 
tics. Examples are (stochastic) attribute- value grammars (Abney, 1997), probabilistic fea- 
ture grammars (Goodman, 1997), head-driven phrase structure grammars (Pollard &: Sag, 
1994), and lexical- functional grammars (Bresnan, 2001). For learning within such frame- 
works, methods from undirected graphical models are used; see the work of Johnson (2003) 
for a description of some recent work. The key difference to LOHMMs is that only nonter- 
minals are replaced with structured, more complex entities. Thus, observation sequences of 
flat symbols and not of atoms are modelled. Goodman's probabilistic feature grammars are 
an exception. They treat terminals and nonterminals as vectors of features. No abstraction 
is made, i.e., the feature vectors are ground instances, and no unification can be employed. 

6. A grammar is in GNF iff all productions are of the form A <— &V where A is a variable, a is exactly one 
terminal and y is a string of none or more variables. 

7. A grammar is in CNF iff every production is of the form A B, C or A ^ a where A, B and C are variables, 
and a is a terminal. 
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Figure 7: (a) Each atom in the logical sequence mkdir(vtlOOx), mv(new*, vtlOOx), 
Is(vtlOOx), cd(vtlOOx) forms a tree. The shaded nodes denote shared labels 
among the trees, (b) The same sequence represented as a single tree. The pred- 
icate con/2 represents the concatenation operator. 



Therefore, the number of parameters that needs to be estimated becomes easily very large, 
data sparsity is a serious problem. Goodman applied smoothing to overcome the problem. 

LOHMMs are generally related to (stochastic) tree automata (see e.g., Car- 
rasco, Oncina, and Calera-Rubio, 2001). Reconsider the Unix command sequence 
mkdir(vtlOOx), mv(new*, VtlOOx), Is(vtlOOx), cd(vtlOOx) . Each atom forms a tree, see 
Figure 7 (a), and, indeed, the whole sequence of atoms also forms a (degenerated) tree, 
see Figure 7 (b). Tree automata process single trees vertically, e.g., bottom-up. A state in 
the automaton is assigned to every node in the tree. The state depends on the node label 
and on the states associated to the siblings of the node. They do not focus on sequential 
domains. In contrast, LOHMMs are intended for learning in sequential domains. They 
process sequences of trees horizontally, i.e., from left to right. Furthermore, unification 
is used to share information between consecutive sequence elements. As Figure 7 (b) 
illustrates, tree automata can only employ this information when allowing higher-order 
transitions, i.e., states depend on their node labels and on the states associated to 
predecessors 1,2,... levels down the tree. 

In the second type of approaches, most attention has been devoted to developing highly 
expressive formalisms, such as e.g. PGUP (Eisele, 1994), PGLP (Riezler, 1998), SLPs (Mug- 
gleton, 1996), PLPs (Ngo & Haddawy, 1997), RBNs (Jaeger, 1997), PRMs (Friedman, 
Getoor, Roller, & Pfeffer, 1999), PRISM (Sato & Kameya, 2001), BLPs (Kersting & De 
Raedt, 2001b, 2001a), and DPRMs (Sanghai, Domingos, & Weld, 2003). LOHMMs can be 
seen as an attempt towards downgrading such highly expressive frameworks. Indeed, apply- 
ing the main idea underlying LOHMMs to non-regular probabilistic grammar, i.e., replacing 
flat symbols with atoms, yields - in principle - stochastic logic programs (Muggleton, 1996). 
As a consequence, LOHMMs represent an interesting position on the expressiveness scale. 
Whereas they retain most of the essential logical features of the more expressive formalisms, 
they seem easier to understand, adapt and learn. This is akin to many contemporary consid- 
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erations in inductive logic programming (Muggleton &: De Raedt, 1994) and multi-relational 
data mining (Dzeroski & Lavrac, 2001). 

8. Conclusions 

Logical hidden Markov models, a new formalism for representing probability distributions 
over sequences of logical atoms, have been introduced and solutions to the three central 
inference problems (evaluation, most likely state sequence and parameter estimation) have 
been provided. Experiments have demonstrated that unification can improve generalization 
accuracy, that the number of parameters of a LOHMM can be an order of magnitude smaller 
than the number of parameters of the corresponding HMM, that the solutions presented 
perform well in practice and also that LOHMMs possess several advantages over traditional 
HMMs for applications involving structured sequences. 
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Appendix A. Proof of Theorem 1 

Let M = (Ti, fi, A,T) be a LOHMM. To show that M specifies a time discrete stochastic 
process, i.e., a sequence of random variables {Xt)t=i,2,..., where the domains of the random 
variable Xt is hb(S), the Herbrand base over S, we define the immediate state operator 
Tjv/-operator and the current emission operator i^Az-operator. 

Definition 4 (Tm- Operator, Em -Operator ) The operators Tm '■ 2'^^^ — )■ 2^*^^ and Em '■ 
2^0^ _^ 2"^"^^ are 

Tm{I) = {Hgb^h \3ip:H ^ B)€ M : BaB G /, HasaH S G^iH)} 

Em{I) = {OoBOHOo \^{p:H ^ B)G M : BaB G /, HaBdc G G^{H) 
and OaBCTHO'o G Gs(0)} 

For each i = 1,2,3,..., the set Ti+^({start}) := rM(Ti/({start})) with 
T^({start}) := TM{{sta.rt}) specifies the state set at clock i which forms a random vari- 
able Yi. The set C/|,j({start}) specifies the possible symbols emitted when transitioning 
from i to i + 1. It forms the variable Ui. Each Yi (resp. Ui) can be extended to a random 
variable Zt (resp. Ui) over hb^: 




otherwise 
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Figure 8: Discrete time stochastic process induced by a LOHMM. The nodes and Ui 
represent random variables over hbs. 



Figure 8 depicts the influence relation among Zi and JJi. Using standard arguments from 
probability theory and noting that 

P{U = U I Z.^i = ..^1, Z. = z.) = ^(^m = ^m,^. = ^. U.) 

z2u, P{Zi+i,ui I Zi) 

and P{Zi+i I Zi) = '^P{Zi+i,Ui | Zi) 

Ui 

where the probability distributions are due to equation (1), it is easy to show that Kol- 
mogorov's extension theorem (see Bauer, 1991; Fristedt and Gray, 1997) holds. Thus, M 
specifies a unique probability distribution over ^^^^^{Zi x Ui) for each t > and in the 
limit t — 7> oo. □ 



Appendix B. Moore Representations of LOHMMs 

For HMMs, Moore representations, i.e., output symbols are associated with states and Mealy 
representations, i.e., output symbols are associated with transitions, are equivalent. In this 
appendix, we will investigate to which extend this also holds for LOHMMs. 

Let L be a Mealy-LOHMM according to definition 3. In the following, we will derive 
the notation of an equivalent LOHMM L' in Moore representation where there are abstract 
transitions and abstract emissions (see below). Each predicate b/n in L is extended to b/n+ 
1 in L'. The domains of the first n arguments are the same as for b/n. The last argument 
will store the observation to be emitted. More precisely, for each abstract transition 

, / N o(vi,...,vi,) 

p : h(wi, . . . , wi) < b(ui, . . . ,Un) 



in L, there is an abstract transition 



p:h(wi,...,Wi,o(vi,...,v(,)) ^b(ui,...,Un,-) 

in L'. The primes in o(v'^, . . . , v(^) denote that we replaced each free ^ variables o(vi, . . . , v^) 
by some distinguished constant symbol, say Due to this, it holds that 

^(h(wi,...,wi)) =/i(h(wi,...,Wi,o(vi,...,v(,))) , (6) 



8. A variable X £ vars(o(vi, . . . is free iff X vars(h(wi, . . . , wi)) U vars(b(ui, . . . , Un)). 
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and L"s output distribution can be specified using abstract emissions which are expressions 
of the form 

1.0 : o(vi, . . . ,Vk) ^ h(wi, . . . ,Wi,o(vi, . . . ,v'^)) . (7) 

The semantics of an abstract transition in L' is that being in some 
state S( G GE'(b(ui, . . . , Uji, -)) the system will make a transition into state 
Sj_|_j € Gs'O^i^i, ■ ■ ■ ,wi,o(v'^, . . . ,v(^))) with probabihty 

P-K^'t+i I li(wi,...,Wi,o(v'i,...,v^)) I o-g/) (8) 

where a^i^ = mgu(S'j, b(ui, . . . ,Uii,_)). Due to Equation (6), Equation (8) can be rewritten 
as 

I h(wi,...,Wl) I CTg/) . 

Due to equation (7), the system will emit the output symbol ot+^ S Gx;'(o(vi, . . . , v^)) in 
state S'(_,_^ with probability 

l-l{0t+i I 0(vi,...,Vk)cTs/^^0-s'J 

where a^/^ ^ = mgu(h(wi, . . . , Wi, o(v'^, . . . , vQ), S'(_|__; ). Due to the construction of L', there 
exists a triple (S^, S(+_? , Ot+i) in L for each triple (S'^, S'j_,__; ), t > 0, in L' (and vise 
versa). Hence, both LOHMMs assign the same overall transition probability. 

L and L' differ only in the way the initialize sequences {{S'^,S'^_^_J ,Ot+i)t=o,2...,T (resp. 
((Sf, S(+i , Ot+i)j=o,2...,T)- Whereas L starts in some state So and makes a transition to Si 
emitting Oi, the Moore-LOHMM L' is supposed to emit a symbol Do in Sg before making a 
transition to S'^^. We compensate for this using the prior distribution. The existence of the 
correct prior distribution for L' can be seen as follows. In L, there are only finitely many 
states reachable at time t = 1, i.e, Pl{qo = S) > holds for only a finite set of ground 
states S. The probability Pl{% = s) can be computed similar to ai(S). We set t = 1 in line 
6, neglecting the condition on Ot-i in line 10, and dropping ^{^t-i \ 0<7BcrH) from line 14. 
Completely listing all states S G 5*1 together with Pl{qo = S), i.e., Pl{(1q = S) : S start , 
constitutes the prior distribution of L' . 

The argumentation basically followed the approach to transform a Mealy machine into 
a Moore machine (see e.g., Hopcroft and Ullman, 1979). Furthermore, the mapping of a 
Moore-LOHMM - as introduced in the present section - into a Mealy-LOHMM is straight- 
forward. 

Appendix C. Proof of Theorem 2 

Let T be a terminal alphabet and N a nonterminal alphabet. A probabilistic context-free 
grammar (PCFG) G consists of a distinguished start symbol S" G iV plus a finite set of 
productions of the form p : X ^ a, where X € N , a € {N U T)* and p G [0, 1]. For all 
X £ N, ^-x-faP = 1- A PCFG defines a stochastic process with sentential forms as states, 
and leftmost rewriting steps as transitions. We denote a single rewriting operation of the 
grammar by a single arrow If as a result of one ore more rewriting operations we are 
able to rewrite /3 G (A^ U T)* as a sequence 7 G {N U T)* of nonterminals and terminals, 
then we write /3 ^* 7. The probability of this rewriting is the product of all probability 
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values associated to productions used in the derivation. We assume G to be consistent, i.e., 
that the sum of all probabilities of derivations S =>* f3 such that f3 T* sum to 1.0. 

We can assume that the PCFG G is in Greibach normal form. This follows from Abney 
et al.'s (1999) Theorem 6 because G is consistent. Thus, every production P G G is of 
the form p : X ^ oYi . . .Yn for some n > 0. In order to encode G as a LOHMM M, we 
introduce (1) for each non-terminal symbol X in G a constant symbol nX and (2) for each 
terminal symbol t in G a constant symbol t. For each production P G G, we include an 
abstract transition of the form p : stack([nYi, . . . ,nYn|S]) ^ stack([nX|S]), if re > 0, and 
p : stack(S) -A stack([nX|S]), if n = 0. Furthermore, we include 1.0 : stack([s]) ^ start 
and 1.0 : end stack([]). It is now straightforward to prove by induction that M and G 
are equivalent. □ 

Appendix D. Logical Hidden Markov Model for Unix Command 
Sequences 

The LOHMMs described below model Unix command sequences triggered by mkdir. To 
this aim, we transformed the original Greenberg data into a sequence of logical atoms over 
com, inkdir(Dir, LastCom), ls(Dir, LastCom), cd(Dir, Dir, LastCom), cp(Dir, Dir, LastCom) 
and mv(Dir, Dir, LastCom). The domain of LastCom was {start, com, mkdir. Is, cd, cp, mv}. 
The domain of Dir consisted of all argument entries for mkdir. Is, cd, cp, mv in the original 
dataset. Switches, pipes, etc. were neglected, and paths were made absolute. This yields 
212 constants in the domain of Dir. All original commands, which were not mkdir. Is, cd, 
cp, or mv, were represented as com. If mkdir did not appear within 10 time steps before a 
command G G {is, cd, cp,mv}, G was represented as com. Overall, this yields more than 
451000 ground states that have to be covered by a Markov model. 

The "unification" LOHMM U basically implements a second order Markov model, i.e., 
the probability of making a transition depends upon the current state and the previous 
state. It has 542 parameters and the following structure: 




mkdir(Dir, com) ^ com. 

end ^ com. 



com ^ 



com. 



Furthermore, for each G G {start, com} there are 




cp(_, _, mkdir) ^ mkdir(Dir,G). 
mv(_, Dir, mkdir) ^ mkdir (Dir, G). 
mv(Dir, _, mkdir) ^ mkdir(Dir,G). 



cd(_, mkdir) ^ mkdir(Dir,G). 
cp(_, Dir, mkdir) ^ mkdir(Dir,G). 
cp(Dir, _, mkdir) <r- mkdir(Dir,G). 



ls(Dir, mkdir) ^ mkdir(Dir,G 



ls(_, mkdir) ^ mkdir(Dir,G 
cd(Dir, mkdir) ^ mkdir (Dir, G 



mv(_, _, mkdir) ^ mkdir(Dir,G). 
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together with for each C G {mkdir, Is, cd, cp, mv} and for each Ci € {cd, Is} (resp. 
C2 e {cp,mv}) 



mkdir(Dir, com) 
mkdir(_, com) 
com 
end 

ls(Dir,Ci) 
ls(-,Ci) 
cd(Dir,Ci) 
cd(_,Ci) 
cp(_,Dir,Ci) 
cp(Dir,_,Ci) 
cp(-,-,Ci) 
mv(_, Dir,Ci) 
mv(Dir, _,Ci) 
mv(_, _,Ci) 



Ci(Dir,C 
Ci(Dir,C 
C7i(Dir,C7 
C7i(Dir,C7 
Ci(Dir,C 
Ci(Dir,C 
Ci(Dir,C 
Ci(Dir,C7 
Ci(Dir,C 
Ci(Dir,C 
Ci(Dir,C 
Ci(Dir,C7 
Ci(Dir,C7 
Ci(Dir,C 



mkdir(_, com) 
com 
end 

ls(From,C2) 
ls(To,C2) 

1S(-,C2) 

cd(From,C2) 
cd(To,C2) 
cd(_,C2) 
cp(From, _,C2) 
cp(_, To,C2) 
cp(_, _,C2) 
mv(From, -^2) 
mv(_, To,C2) 
mv(_, _,C2) 



C2(From, 
C2(From, 
C2(From, 
C2(From, 
C2(From, 
C2(From, 
C2(From, 
C2(From, 
C2(From, 
C2(From, 
C2(From, 
C2(From, 
C2(From, 
C2(From, 
C2(From 



To,C). 
To,C). 
To,C). 
To,C). 
To,C). 
To,C). 
To,C). 
To,C). 
To,C7). 
To,C). 
To,C). 
To,C). 
To,C). 
To,C). 
To,C). 



Because all states are fully observable, we omitted the output symbols associated with 
clauses, and, for the sake of simplicity, we omitted associated probability values. 

The "no unification" LOHMM N is the variant of U where no variables were shared 
such as 



mkdir(_, com) ^ cp(From, To,C). 

com ^ cp(From, To,C). 
end ^ cp(From, To,C). 



ls(_, cp) ^ cp(From, To,C). 

cd(_, cp) ^ cp(From, To,C). 

cp(_, _, cp) cp(From, To,C). 

mv(_, _, cp) cp(From, To,C). 



Because only transitions are affected, has 164 parameters less than [/, i.e., 378. 



Appendix E. Tree-based LOHMM for mRNA Sequences 

The LOHMM processes the nodes of mRNA trees in in-order. The structure of the LOHMM 
is shown at the end of the section. There are copies of the shaded parts. Terms are 
abbreviated using their starting alphanumerical; tr stands for tree, he for helical, si for 
single, nuc for nucleotide, and nuc_p for nucleotide_pair. 

The domain of ^Children covers the maximal branching factor found in the data, i.e., 
{[c], [c, c], . . . , [c, c, c, c, c, c, c, c, c]}; the domain of Type consists of all types occurring in 
the data, i.e., {stem, single, bulgeS, bulgeS, hairpin}; and for Size, the domain covers 
the maximal length of a secondary structure element in the data, i.e., the longest sequence 
of consecutive bases respectively base pairs constituting a secondary structure element. 
The length was encoded as {n^(0), n^(0), . . . , n^^(O)} where n"^(0) denotes the recursive 
application of the functor n m times. For Base and BasePair, the domains were the 4 bases 
respectively the 16 base pairs. In total, there are 491 parameters. 
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Copies for trCid, [c] , [Pa - [C]|R]), tr(Id, [c, c], [Pa - [C]|R]), Copies for tr(ld, [c] , [Pa - [Cl, C2|Cs]|R]), tr(ld, [c, c], [Pa - [Cl, C2|Cs]|R] | 
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